The QCD collisional energy loss revised 
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It is shown that to leading order the QCD collisional energy loss reads dE/dx ~ OL(m 2 D )T 2 . Com- 
pared to prevalent expressions, dE B /dx ~ a 2 T 2 \n(ET/m 2 D ), which could be considered adaptions 
of the (QED) Bethe-Bloch formula, the rectified result takes into account the running coupling, 
as dictated by renormalization. As one significant consequence, due to asymptotic freedom, the 
collisional energy loss becomes independent of the jet energy E. Some implications with regard to 
heavy ion collisions are pointed out. 
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One of the key arguments for the creation of a 'new 
state of matter' in heavy ion collisions at RHIC is the 
observed jet quenching [I], which inter alia probes the 
jet's energy loss in the traversed matter. In a quark gluon 
plasma (QGP) there are two effects causing a jet to lose 
energy: elastic collisions with deconfined partons [2, 3], 
or induced gluon radiation [4, 5, 6]. Presuming a dom- 
inance of the second mechanism, experimental findings 
have often been interpreted in terms of a purely radia- 
tive loss. However, the data-adjusted parameters (either 
q or dNg/dy, depending on the approach) are found to be 
considerably larger than theoretically expected or even 
in conflict with a strong constraint from dS/dy (see e. g. 
[7]) - which calls for a collisional component of the en- 
ergy loss. The effect of collisions (as estimated within 
the framework [2, 3]) might actually be larger than con- 
ceded for a long time [8, 9]. The fact that such estimates 
depend crucially on the assumed value of the coupling 
should motivate us to scrutinize the principal question 
' What is a?'. Aside from its phenomcnological relevance 
it will also lead to interesting theoretical insight. 

Following Bjorkcn's intuitive considerations [2], con- 
sider the propagation of a jet through a static QGP at 
a temperature T A, where the coupling is small. Its 
mean energy loss per length can be calculated from the 
rate of binary collisions with partons of the medium, as 
determined by the flux and the cross section, 
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Here p s = d s n s is the density of scatterers, with d g = 16 
and d q = Ylnj being the gluon and quark degeneracies 
for rif light flavors, and n±(k) = [exp(fc/T) ± 1] 1 in 
the ideal gas approximation. Furthermore, $ denotes 
a dimcnsionless flux factor, t the 4-momentum transfer 
squared, and u> = E — E' the energy difference of the 
incoming and outgoing jet. Focusing on the dominat- 
ing scatterings with small momentum transfer, the cross 
sections read 
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with C a 



|, C qg = 1, and C gg = §. For E and E' 
much larger than the typical momentum of the thermal 



scatterers, k ~ T, the relation of t to the angle 9 between 
the jet and the scatterer simplifies to 

i = -2(l-cos0)fej, (3) 

and the flux factor can be approximated by $ = 1 — cos 8. 
At this point, Bjorken integrated in Eq. (1) 
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imposing both an IR and UV regularization. The soft 
cut-off is related to the Debye mass, ^1 = p 2 ~ 
m D ~ aT 2 , describing the screening of the exchanged 
gluon in the medium. The upper bound of |t| was 
reasoned to be given by the maximum energy trans- 
fer: very hard transfers, say uj E, effectively do 
not contribute to the energy loss; in this case the en- 
ergy is collinearly relocated to the scatterer. Assum- 
ing u max = E/2 implies t\ = —(1 — cos0)kE, hence 
dEf/dx = T7a 2 Y, s QsJ k sk- 1 p s \n((l~cose)kE/fi 2 ). 
Replacing now, somewhat pragmatically, the logarithm 
by ln(2(k)E/(j> 2 ) and setting (k) — ► 2T, Bjorken obtained 
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which differs for quark and gluon jets only by the pref- 
actor. This expression can be regarded as a rclativistic 
adaption of the (QED) Bcthc-Bloch formula [10], which 
describes the ionization/excitation energy loss of charged 
particles in matter as determined by the scatterer den- 
sity, and with the logarithm reflecting the kinematics and 
the long-range Coulomb-type interaction. 

There are various refinements of Bjorken's 'practical' 
calculation, aiming at the precise determination of the 
cut-offs. Worth accentuating is the approach of Braaten 
and Thoma [3] who studied, within HTL perturbation 
theory, the propagation of a fcrmion through a QED 
plasma (and applied their method also to QCD). For light 
quarks, the result reproduces the generic form (5), with 
p — * mjj in the logarithm and the factor 4 replaced by 
some function of [18]. 

A remark concerning a pragmatic usage of such 
'Bjorkcn-type' formulae seems apposite here. Applied to 
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experimentally relevant temperatures and rather low E, 
a formally resulting negative loss has been interpreted, at 
times, as an energy transfer of the thermal medium to the 
'soft jet'. This interpretation, however, is untenable in 
the given framework: the jet always loses energy in a col- 
lision, u) > 0, cf. (3). A negative result for dE B jdx is de 
facto the consequence of interchanged boundaries in the 
integration (4). Since /i 2 is the minimal \t\ by definition, 
dE B jdx should instead be set to zero for AE < fi 2 /T. 
This concern, though, will prove irrelevant by the follow- 
ing considerations. 

In Bjorken's derivation, a is a fixed parameter. Con- 
ceptionally, such a tree-level approximation may be ap- 
propriate for QED. The strong interaction, however, is 
known to vary considerably over the range of scales 
probed, e.g., in heavy ion collisions. Thus, in QCD one 
should study quantum corrections to the tree level pro- 
cesses, whose renormalization will specify unambiguously 
the value of 'the' coupling. 

For the sake of transparency of the argument, consider 
first the analogous case of electron scattering in mass- 
less QED. There are three types of (divergent) loop cor- 
rections to the ^-channel tree-level process, see e.g. [11]. 
First, the exchanged photon is dressed by a self-energy. 
Then, encoded in the quantity Zi, there arc vertex cor- 
rections and finally, via the field strength renormalization 
Z2, self-energy corrections to the external fermions. Yet, 
due to the identity Z\ = Z2, in QED the renormalized 
coupling is determined only by the boson self-energy. 

It is appropriate to renormalize the theory (i.e, to 
fix its parameters) by a scattering experiment at T = 
0. The relevant part of the matrix element leading 
to the (vacuum) cross section corresponding to (2) is 
a/(P 2 — II vac (P 2 )). Here a denotes the bare coupling, 
and II vac is the unrenormalized boson self-energy at 
T = 0. In dimensional regularization, II vac (P 2 ) = 
ab [e _1 - \n(-P 2 /n 2 )] P 2 , where b = 4ir(3 and (3 is 
the leading coefficient of the /3-function. For a specific 
P 2 = t r , the matrix element reads explicitly 



1 a ce(t r ) 

U 1 - ab [e- 1 - ln(-i r /^ 2 )] = t r ' 

A measurement then specifies the renormalized coupling 
a(t r ) at the scale t r (as introduced on the rhs) which is 
related to the (infinite) bare coupling by 

a' 1 (U) = a- 1 - b Q [e- 1 - ln(-t r /^ 2 )] . (6) 

An equivalent relation holds for the coupling a(t) at 
an arbitrary scale t, consequently a~ l (t) = a _1 (t r ) + 
6q \n(t/t r ) or, in a common alternative form, 



a(t) = [b ln(|t|/A 2 )] 
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It is underlined that the momentum dependence of the 
renormalized ('running') coupling is fully specified by its 
value at a certain scale t r or, equivalcntly, by the param- 
eter A. 



In a (thermal) medium, the boson self-energy has the 
generic structure 
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ab [(e- 1 - H-P 2 /^ 2 )) P 2 + r(p ,p)} 



where the finite 'matter' contributions /' differentiate 
transverse and longitudinal modes (i =t,l). Then, utiliz- 
ing (6), the in- medium scattering matrix can be rewritten 
in terms of the renormalized coupling, 
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This distinct form, where the divergent vacuum contribu- 
tion of the self-energy is 'absorbed' in the running cou- 
pling, elucidates that the matrix element depends only 
on the physical parameter A and the matter part of the 
self-energy, IP mat = a(P 2 )b f l ■ Thus the effective IR 
cut-off for the energy loss is, as expected, related to the 
Debye mass [19]. The main emphasis here, though, is on 
the renormalized coupling in Eq. (8) and, consequently, 
also in the resulting differential cross section: the scale 
of the running coupling is set by the virtuality P" = t. 

It is physically intuitive that this fact is generic. Thus, 
instead of a detailed analysis of loop corrections in QCD 
(which is more complex), it is useful to invoke a more 
comprehensive argument. The vacuum differential cross 
section, as a quantity with an unambiguous normaliza- 
tion, obeys a Callan-Symanzik equation with 7 = [11], 
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where g = (47ra(M 2 )) is related to the coupling at a 
given renormalization point M. The general solution of 
this linear partial differential equation is a function h(x), 
whose argument x = g(M) satisfies Mdg/dM + (3(g) = 0. 
To investigate, in this line of thoughts, the dependence on 
the momentum scale Q, with Q 2 = q 2 = t, note that, in 
the limit of small t, the most general form of the cross sec- 
tion is da/dt = S(Q/M,g)/t 2 . Consequently, the func- 
tion S obeys [—Qdq + (3(g) d g ]S = (mind the minus 
sign). The corresponding characteristic equation, with 
0(d) = ~flog 3 at leading order, then leads readily to the 
running coupling as introduced above. In other words, 
renormalization group invariance implies that loop cor- 
rections to the differential cross section can be 'absorbed' 
in the tree level expression by replacing a — > a(t). 

A running coupling, as given by Eq. (7) with 60 = 
(11 - §n/)/(47r) in QCD, alters the integral (4) [20], 
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Opposed to the expression (4), the weighted cross section 
is UV finite - due to the asymptotic weakening of the 
strong interaction, large-t contributions are suppressed. 
For hard jets, the integral (10) becomes independent of 
the energy E, i. e., it is then controlled solely by the cou- 
pling at the screening scale. The necessary condition 
| ti| ~ TE p? ~ m 2 D ~ aT 2 , is, for weak coupling, 
actually much less restrictive than the previous assump- 
tion E 3> T to simplify the kinematics. Anticipating an 
extrapolation of the final result for dE/dx to larger cou- 
pling, where the hierarchy of scales y/aT <C T breaks 
down, it is mentioned that mjj ^ 3T (see, e.g., [12]). 
Hence from this perspective, the collisional loss could 
become ^-independent already for jet energies exceeding 
several GeV's - provided, of course, that the perturbative 
framework gives at least a scmi-quantitativc guidance at 
larger coupling, which will be advocated below. 
The jet's collisional energy loss approaches, thus, 
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In the ideal gas limit, the remaining integration yields 
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which departs radically from previous expressions as 
Bjorken's (5). Aside from the modified cut-off depen- 
dence discussed above, the QCD collisional loss is pro- 
portional to a (instead of a 2 ). It is highlighted that 
the considerations also show that the relevant scale for 
the coupling is the (perturbatively soft) screening mass 
rather than a 'characteristic' thermal (hard) scale ~ T, 
as commonly presumed. 

In order to quantitatively compare Eq. (12) to previous 
estimates it is necessary to specify parameters, namely A 
in (7) and the cut-off fj,, i.e., the Debye mass. Similar 
renormalization arguments as above lead in [12] to an 
implicit equation for mo', to leading order 
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whose solution can be given in terms of Lambert's func- 
tion. Obviously, this mended perturbative formula is jus- 
tified strictly only at temperatures T ^> A. Notwith- 
standing this, it is found in quantitative agreement with 
lattice QCD calculations down to T w 1.2T C [21]. 

It may come as a surprise that a leading order for- 
mula reproduces non-perturbative results [22]. Thus, 
it is worth emphasizing that the adjusted parameter 
A = 205 MeV for rif = 2 is right in the expected ball- 
park, enervating a possibility of an uninterpretable 'over- 
strained' fit. Moreover, the 1-loop running coupling (7) 
with the same A reproduces lattice calculations for an- 
other quantity, namely the QCD vacuum potential V{r), 
in fact up to large distances corresponding to a(r~ 2 ) m 1 
[12]. Although at first sight rather different, the poten- 
tial at T = and the Dcbyc mass in the medium arc 



closely related by the renormalized t-channel scattering 
discussed above - as is the collisional energy loss. In 
other words (tidying the order of the arguments): one 
can renormalize the theory at T = (i.e., determine 
once and for all A from V(r)), verify the applicability 
of the perturbative approach for larger couplings as rele- 
vant near T c by successfully calculating mo (T) , and then 
make predictions for dE/dx. 

With this justification, an extrapolation of Eq. (12) 
as presented in Fig. 1 might be not too unreasonable. 
Assumed here is T c = 160 MeV and p 2 = [\,2]m 2 D to 
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FIG. 1: Light quark collisional energy loss: Eq. (12) vs. the 
prevalent expression (5) (which yields negative values for very 
large T) for representative jet energies. For details see text. 

estimate the uncertainty due to the IR cut-off. Shown for 
comparison are results from (5); here fi = mo (likewise 
from (13)) albeit a in the prefactor (unjustified, but as 
often presumed) fixed at the scale Qt = 2irT. It turns 
out that already near T c , the estimates from formula (12) 
exceed those from (5), even for rather large values of E. 

For phcnomcnological implications it is instructive to 
take further into account a main effect of the strong in- 
teraction near T c . From the distinct decrease of the QGP 
entropy seen in lattice QCD calculations [14], one can, on 
general grounds, infer a similar behavior for the number 
densities p s . In the framework of the quasiparticlc model 
[15], the ideal distribution functions in Eq. (11) are to 
be replaced by d s n±{\J m 2 (T) + fc 2 ), where the effective 
coupling in the quasiparticle masses, m 2 cx a c s(T)T 2 , 
is adjusted to lattice results for the entropy. As shown 
in Fig. 2, near T c - despite the strong interaction - the 
plasma becomes transparent due to the reduced number 
of 'active' degrees of freedom. Qualitatively, this charac- 
teristic behavior is in line with a rather sudden transition 
from small to large jet quenching when going from SPS 
to RHIC energies. 

In conclusion, it has been demonstrated in the context 
of thermal field theory that renormalization does not only 
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FIG. 2: Influence of reduced scatterer density near T c on the 
gluon energy loss. Below T w 1.2T C , Eq. (13) over-estimates 
the Debye mass as obtained in lattice QCD; hence dE/dx 
could be slightly larger than depicted by the dotted line. 

dictate the value of the coupling for a given quantity, but 
that the running of the coupling can also influence the 
parametric behavior of results. For the QCD collisional 
energy loss, the relevant scale in a(t) is the (perturba- 
tively soft) screening mass mrj ~ \faT (instead of T, as 
commonly presumed). The increasing coupling at soft 



momenta leads to a parametric enhancement compared 
to previous calculations, see Eqs. (12) vs. (5). On the 
other hand, due to asymptotic freedom, the collisional 
energy loss becomes independent of the jet energy E. 
Thus, the asymptotic behavior of the underlying interac- 
tion makes the energy loss qualitatively different in QCD 
and QED (where an analog of Eq. (5) indeed holds). 

Except very near T c , Eq. (12) suggests a larger colli- 
sional energy loss than previously estimated [2, 3], This 
finding can be interpreted as a facet of the 'strongly 
coupled' QGP (sQGP), which is characterized by large 
cross sections. In fact, a = Jdtda/dt with running 
coupling can actually be an order of magnitude larger 
than expected from the widely used expression cr a B x ex 
cv 2 (Qy)//i 2 . Thus, the present approach gives a consis- 
tent and simple explanation of phcnomenologically in- 
ferred er ~ O(10)mb [16, 17]. 

Close to T c , though, the particle density is known to 
be substantially reduced. This implies that, irrespective 
of the strong coupling, the sQGP becomes transparent 
near the (phase) transition. Such a distinct temperature 
dependence of the energy loss should be observable. The 
quantification of this effect (including a discussion impli- 
cations of the running coupling for the radiative energy 
loss) will be the subject of a forthcoming study. 
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